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" Abstract 

■ Spectrum of low- lying five-quark configurations with strangeness quantum number S = — 3 

7— i ; 

qq ■ and negative parity is studied in three kinds of constituent quark models, namely the one gluon 

o: 

£S) , exchange, Goldstone Boson exchange, and instanton-induced hyperfine interaction models, respec- 

tively. Our numerical results show that the lowest energy states in all the three employed models 
^ ■ are lying at ~1800 MeV, about 200 MeV lower than predictions of various quenched three-quark 

h; 

models. In addition, it is very interesting that the state with the lowest energy in one gluon 
exchange model is with spin 3/2, but 1/2 in the other two models. 
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I. INTRODUCTION 



Up to now, the experimental information on the f2 hyperon excitations is still very 
poor [l]. There are only four f2 states discovered {2^5], namely the ground state f2(1672), 
and three excitations ^(2250), ^(2380) and fi(2470). Among these four states, only ^(1672) 
has been clarified to have spin parity quantum number | , while quantum numbers for the 
other three have not been pinned down yet. 

Although there has not been any further experimental evidence about Q* since 1990s, 
theorists are always interested in spectrum of Q excitations, which has been investigate^ . 
within various approaches, such as constituent quark models based on one gloun exchange 



Q], Goldstone Boson exchange 
large 1/N C expansion of QCD 
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9j and instanton-induced 



10] interactions, respectively, 



iMlq. bound state approach of Skyrme model 16j, meson- 
baryon chiral dynamics 17| . etc. Most of the theorists have tried to compare their predictions 
with data, while information from experiments is so poor that definite comparison seems 
to be impossible. However, to give more guidance for future experimental measurements, 
systematical studies of Q* spectrum in different theoretical approaches are valuable. 

Recently, the traditional constituent quark model is developed to include higher Fock 
component s,it is shown that there should be notable five-quark components in the baryon 



resonances 



18H27J. So it immediately occurs to us that the five-quark components in fi* 



should be very important, because all the three valence quarks in Q* are strange ones, then 
the five-quark configurations with qq pairs (we denote QQ = qq and QQ = ss for light and 
strange quark- ant iquark pairs, respectively.) must play very special roles in the properties of 
Q*. Since the qq pair has different flavor with valence strange quarks, it has two advantages 
for the study of five-quark components. First, there is not any Pauli blocking effect for the q, 
which would result in lower excitation energy for the new excitation mechanism of Q states 
by pulling out a qq pair and make the five-quark components larger in Q* than in N*, A*, 
etc. Secondly it simplifies the model calculation of the five quark system. 

To specify the roles of s 3 QQ in Q*, as the first step, we investigate the spectrum of low 
lying s 3 QQ configurations with negative parity in this manuscript, employing the constituent 
quark model complemented by three different kinds of hyperfine interactions between quarks, 
i.e. the one gluon exchange (OGE), Goldstone Boson exchange (GBE) and instanton-induced 
(INS) force interactions, respectively. In order to get more accurate information about the 
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five-quark configurations with light qq pairs, the S77(3) (&SU(2) Cg>0(3) breaking corrections 
are taken into account in our model. 

The present manuscript is organized as follows. In Section [III we present our theoretical 
framework, which includes explicit forms of the employed hyperfine interactions between 
quarks. Numerical results for spectrum of the studied five-quark configurations in our model 
are shown in Section IIHI Finally, Section IIVI contains a brief conclusion. 



II. THEORETICAL FRAMEWORK 

In the constituent quark model, Hamiltonian for a five-quark system with a QQ pair can 
be expressed as: 

5 

H = Ho + Hhyp + Y.mi, (1) 
i=i 

where denote the constituent masses of the quarks, Hh yp is the hyperfine interaction 
between quarks, which is often treated as perturbation, and H Q the Hamiltonian concerning 
orbital motions of the quarks which should contain two parts, namely the kinetic term and 
confining potential of the quarks. Both of the orbital Hamiltonian H Q and the hyperfine 
interaction Hh yp for three-quark system have been intensively discussed in literatures, here 
we will develop them to five-quark system explicitly in Sec. Ill Al and III B[ 



A. The orbital Hamiltonian for a five-quark system 
In this manuscript, the employed form of H Q is as follow 

where P cm and M denote the total momentum and mass of the five-quark system. The 
confining potential V con f(rij) is taken to be the commonly accepted color confinement 

Vconfir^) = -hf ■ Xf [C{n - r 3 f + Vb] • (3) 

Here A&s are the Gell-Mann matrices in the 577(3) color space, C denotes a confinement 
strength constant, and Vq represents the unharmonic part of V con f which is treated as a 
constant in this manuscript. Explicit calculations show that the matrix elements of Xf ■ Xf 



between both of QQ and QQ pairs in a five-quark system are the same value —4/3, which 
is half of the matrix element between a QQ pair in a three-quark system. Therefore, in the 
case of the QQ in a five-quark system being ss pair, H Q reduces to the following form 

H o = J2 (ft/2m s + + 5V Q , (4) 

where £i represent the Jacobi coordinates for a system of 5 constituents which are defined 
as follow 




Si 

\ + I 2 
1 5 

Rem = T r i i (5) 
5 i=l 

and Pj the corresponding momentums. Then the eigenvalue of if D should be 

E = (N + 6 )w + 5F , (6) 



where u denotes the oscillator frequency which is defined as cu = y5C/m s , and N = 2N r + L 
is the number of excited quanta of the oscillators. In this manuscript, we consider the s 3 QQ 
configurations with N r = L = 0, namely all the quarks and antiquarks in the five-quark 
configurations are in their ground states, then the energy (jHJ) reduces to E Q = 6u + 5Vq. lo 
can be determined in the following way, as we know, the oscillator frequency U3 for a three 



light quarks system is = jQC/m with m the constituent mass of the light quark, then 



we can obtain that cu = \/ ^-^3- And U3 can be inferred from the empirical radius of the 



proton which leads to ^muj 3 ~ 246 MeV, or mass splitting between nucleon and iV(1440) 
which yields oj% = 251 MeV. We take the latter value in this manuscript. 

Eigenfunction of H Q in the case of iV = is just a completely symmetric wave function 
regarding the four Jaccobi coordinates, then we can construct the complete wave functions 
of s 3 QQ configurations in the following way: first, color wave functions of the four-quark 
subsystems in any QQQQQ systems must be [211]c, to combine with the anti-quark to 
form a color singlet, consequently, the flavor-spin symmetry of the four-quark subsystems 



: a 
3, 



in the s 3 QQ configurations can only be [3 1]fs- Moreover, the flavor wave functions o 
sssq system can be completely symmetric [4]jr or mixed symmetric [31}f- As shown in 
four kinds of four-quark wave functions can meet the above requirements. Coupling these 
four-quark systems to the antiquark, one can obtain two five-quark configurations with ss, 



and seven configurations with qq. The total spin for these configurations can be 1/2, 3/2 
and 5/2, as shown in Tabled] Explicit wave functions for these configurations can be easily 



derived from the ones given in 
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Finally, in the case of QQ being a light quark- ant iquark pair, we have to take into account 
the SU(3) breaking corrections. As shown in 29J, we can use the Hamiltonian (JJ|) as an 
approximation which can be complemented by a flavor dependent perturbative term H 



h: 



4 / — jO . — jO 



E 

8=1 



P- 



2m 



5(3m s + 2m) 



5 



iq ■ 



(7) 



here m is the constituent mass of the light quark, 5i q is a flavor dependent operator with 

eigenvalue 1 for light quark and for strange quark. In the present case, all of the quarks 

and antiquark are in their ground states, one can get that the matrix elements of H between 

all of the studied configurations are the same one 

6 ,m s 
5 m 



(K) 



(8) 



B. The hyperfine interaction between quarks 



As stated in Sec III there are nine different configurations studied in this manuscript. 
If the approximation that H^ yp = is applied, the five-quark configurations with qq and ss 
should be two categories of degenerate states. To calculate the mass splittings of the degen- 
erate configurations, explicit perturbative hyperfine interactions are needed. Here we em ploy 
three different kinds of spimspjn forces, which are mediated by one gluon exchange 30|-|33], 
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34 



37[, respectively, 



Goldstone Boson exchange Is |9f and instanton-induced interactions 
as hyperfine interactions between quarks. These interactions between both of QQ and QQ 
pairs have been explicitly given in literatures, here we just present a very brief review, and 
apply them to a five-quark system. 

In the OGE model, hyperfine interaction between quarks can be expressed as chromo- 
magnetic form, which has been intensively used in the studies of multiquark configurations 



38 



-|41|. Here we employ the form given in 



40 



4l| 



ttOGE 

hyp 



(9) 



where c?j are the Pauli spin matrices, Af the Gell-Mann SU(3)c color matrices, and C^j 
a flavor dependent chromomagnetic interaction strength operator. If the subscript i or j 



indicates an antiquark, the following replacement should be applied: \ c — > —\ c *. Once the 
overall symmetry is taken into account, the hyperfine interaction ([9]) should reduce to 

= -6Ci l2 Af • ■ B 2 + 4C 4 , 5 A? • If a 4 ■ a 5 (10) 

The resulting matrix elements of the color, spin and flavor dependent chromomagnetic 
strength operators are given explicitly in Appe ndix lAl And values of interaction strength 
constants Cqq and Cqq can be found in 401442 1 . 



In Ref. [9[, the GBE model is applied to calculate spectrum of the five-quark configura- 
tions with quantum numbers of N, A, E, H and Q resonances. While in that paper, the 
hyperfine interaction is of a SU(3) symmetric form, namely, the coupling strength of n, K 
and rj meson exchanges are taken to be a same one. In present case, we take the SU(3) 



broken form of which is given in 

4 

H hy B p E = - Yl C lj^i ' ^j^i ' °3i 

where Af are the Gell-Mann matrices in flavor space, and Cf^ a flavor dependent operator 
for strength of a meson M exchange between the ith and jth quarks. Notice that in the GBE 
model 8|, hyperfine interaction between quark and antiquark is assumed to be automatically 
included in the GBE interaction, so the spin-spin interaction in Eq. ( TTTj) is restricted 

to the four-quark subsystem. After taking into account the overall symmetry, reduces 
to the following form 

Hgf = -6^ 2 Af ■ Xfo ■ a 2 , (12) 

the matrix elements of the spin operator o\ ■ a 2 are the same as those in the OGE model 
which are given in Appendix and we list the matrix elements of the flavor dependent 
operator Cf 4 2 Xf ■ X F in Appendix [B] Notice that we have to treat properly the three subsets 
of the r) meson exchange between two quarks, namely, the exchange interaction for pairs of 
light qnaiks(C uu and C dd ), and two strange quarks (C ss ). It's proposed by Glozman and 
Riska in that the strengths of the meson exchange interactions may have the following 
relations 

_ rn try* 

C uu = (jdd = C n. C K = git. C ss = / , (13) 



and the phenomenological value for C n is around 20~30 MeV 



BE 
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Finally, the instanton-induced hyperfine interaction between quarks is of the form 35M37| 



(14) 



where the first term is the operator for interaction between two quarks, and the second for 
quark-antiquark interaction, pW are the projectors on flavor antisymmetrical QQ' states, 
ps=x p ro j ec tors on spin x states, and P c,y the projectors on color representation of 
dimension y, respectively, g and g' are the strength constants of interactions between two 
light quarks, and one light quark and one strange quark, respectively. Finally, is a flavor 
projector operator, its eigenvalue is -g for qq and -g for sq. Notice that g also couples isospin 
QQ states. Matrix elements of spin and color operators are given in Appendix The 
values of g and g' can be found in Ref j^j], notice that a factor 6 is included in g and g' 
in [42]. 



III. NUMERICAL RESULTS 

As stated in Sec [ITJ, parameters in the present manuscript are the constituent masses of 
the light and strange quarks m and m s , the unharmonic part of the confinement potential 
Vq, the oscillator parameter u, and the strength constants in the employed three kinds of 
hyperfine interactions. We just take the parameters for hyperfine interactions to be the ones 
given in the literatures listed in Sec UH and the constituent masses of quarks and Vq to be 



the empirical ones given in |42j, and the oscillator parameter u is then ~ 196 MeV. All the 
values for parameters are listed in Table [TTl 

The procedure to obtain numerical results in our model can be via three steps: at first, 
one has to calculate the eigenvalues of the orbital Hamiltonian H Q , which can be obtained 
by Eqs. ([6]) and (151); secondly, one should calculate the matrix elements of the perturbative 
hyperfine interactions, as we can see in Table [Ij matrices for the spin 1/2, 3/2 and 5/2 
configurations are with 4, 4 and 1 dimensions, respectively; finally, energies of the studied 
states should be produced by diagonalization of the obtained matrices in the second step. 

In a given model, all of the studied five-quark configurations should have a same energy 
if we don't consider the hyperfine interactions between quarks and difference between con- 
stituent masses of light and strange quarks. Since the unharmonic parameter Vq is taken 
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to be different in the three hyperfine interaction model 42j, the obtained numerical values 
Eq = Y%=i m i + E for the five-quark configurations with qq in the OGE, GBE and INS 
models are 2196, 1891 and 2171 MeV, respectively, and those for the configurations with ss 
are 2436, 2131 and 2411 MeV in the three models, respectively. 

In the calculations of elements of the hyperfine interactions matrices, we label the con- 
figurations listed in Table [J as follow 

1 : S 3 g([4] x [211] c [31] F5 [31] F [22] 5 )®g, 

2 : S 3 g([4] x [211] c [31] F5 [31] F [31] 5 )®g, 

3 : S 3 g([4] x [211] c [31] F5 [4] F [31] s )®g, 

4 : S 4 ([4] x [211] c [31] F5 [4] F [31] 5 )®s, (15) 
for the configurations with total spin 1/2, and 

1 : S 3 g([4] x [211] c [31] F5 [31] F [31] 5 )®g, 

2 : S 3 g([4] x [211] c [31] F5 [31] F [4] 5 )®g, 

3 : S 3 g([4] x [211] c [31] F5 [4] F [31] s )®g, 

4 : S 4 ([4] x [211] c [31] FS [4] F [31] s )(8)s, (16) 

for the configurations with total spin 3/2. We give the numerical results for the energies of 
these configurations in the next two subsections. 



A. Numerical results without corrections of configuration mixing 



Explicit calculations of the matrix elements of hyperfine interactions between quarks 
in the three models lead to the following matrices for the spin 1/2 and 3/2 five-quark 
configurations 

/ 2235.0 -139.6 10.8 0.0 \ /2223.4 328.9 6.6 0.0 \ 



cOGE 
C l/2 



-149.6 2365.4 -25.6 0.0 
10.8 -25.6 2373.7 0.0 
V 0.0 0.0 0.0 2654.7/ 



cOGE 
' °3/2 



328.9 2095.0 -68.0 0.0 
6.6 -68.0 2333.7 0.0 
V 0.0 0.0 0.0 2517.1 ) 
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GBE 
/2 



cINS 
C l/2 



( 1833.6 0.0 0.0 0.0 \ 

0.0 1896.6 -16.2 0.0 

0.0 -16.2 2010.0 0.0 

V 0.0 0.0 0.0 2161.6/ 

/ 1928.0 -121.5 -30.4 -33.3 \ 

121.5 1908.8 30.4 33.3 



cGBE 
> c 3/2 



/ 1896.6 0.0 
0.0 1990.2 



-16.2 



0.0 



-16.2 

0.0 
2010.0 



0.0 \ 
0.0 
0.0 



cINS 
) °3/2 



\ 0.0 0.0 0.0 2161.6/ 
/ 2052.0 -113.3 -60.7 -66.7 \ 

-113.3 2250.0 191.9 210i 



(17) 



-30.4 30.4 2230.3 47.1 -60.7 191.9 2159.0 188.5 

V -33.3 33.3 47.1 2411.0/ V -66.7 210.8 188.5 2411.0/ 

The results are in unit of MeV. Notice that values of Eq are already included in the diagonal 
terms of the matrices. For the only configuration with spin 5/2, the numerical results in 
OGE, GBE and INS models are 2492, 1990 and 1987 MeV, respectively. 

As we can see in Eq. (IT7|) . in the OGE model, it's very interesting that the diagonal 
matrix element with lowest energy in is about 140 MeV higher than that in £$% E . 

As we know, the diagonal matrix elements £^ yp (i,i) in the above matrices just represent 
the energies of the ith five-quark configuration as labeled by Eqs. f fT5|) and ( Tl6]) . so we 
can conclude that the lowest energy state in OGE model is with spin 3/2, if mixing of 
the configurations with same spin are not taken into account. It's because of that the 
contributions from QQ hyperfine interactions to spin 1/2 states are larger than those to 3/2 
states. This is the most important difference between the OGE and the other two models. 
In addition, as we can see in ffl7|) . some of the non-diagonal terms Mg yp (i,j) are very large, 
which should lead to strong mixing between different configurations. While note that the 
configuration with ss pair doesn't mix with the other three configurations in OGE model. 

In the GBE model, there are only two nonzero non-diagonal matrix elements in both of 
£y BE and £y BE ' ■ These nonzero elements are caused by the non- vanishing matrix elements 
of the flavor dependent operator between flavor wave functions [31] f and [4]^. But the 
nonzero £g BE (i,j) are much smaller than the diagonal matrix elements £g BE (i,i), so the 
mixing of the configurations with same spin should be very weak in GBE model. One may 
notice that an obvious difference between OGE and GBE models is that QQ interactions 
are assumed to be included in meson-exchange between quarks in GBE model, but in OGE 
model, QQ interactions which lead to large non-diagonal matrix elements are independent 
with hyperfine interactions between quark pairs. 

In the INS model, the most significant feature is that the five-quark configurations with 
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ss can mix with those with light qq, as we can see in Eq. ( TTTj) . all the non-diagonal ma- 
trix elements S^ NS (i,4:) with i = 1,2,3 are not 0, this is caused by the nonzero matrix 
elements (qq\gi,5\ss) . Moreover, the matrix element (uu\gi^\dd) is also not 0. In the present 
manuscript, the studied five-quark configurations are with quantum numbers of Q*, whose 
isospin should be 0, consequently, in the case of QQ = qq, wave function should be the 
isocalar one 1/ ' \/2(\sssuu) + \sssdd)). Therefore, the QQ interaction between u and d com- 
ponents should also contribute to the total energy. 



B. Numerical results with configurations mixings 

Diagonalization of matrices (fl7|) leads to the energies for the studied five-quark configu- 
rations as shown in Table ITO1 compared to what obtained by Helminen and Riska js] using 
SU(3) symmetric version of GBE model, and the wave functions for spin 1/2 and 3/2 states 
correspond to the energies in Table HTT1 are of course linear combinations of the configurations 
listed in Table [U and corresponding coefficients for the combinations are given in Table IIVI 
and|Vl 

As we can see in Table IHIt in OGE model, most of the obtained energies after diago- 
nalization of hyperfine interaction matrices are lower than the diagonal matrix elements in 
Eq. f|T7|) . it means that mixing of the configurations decreases the energies in OGE model. 
While in any case, spin of the state with lowest energy is still 3/2, which is about 328 
MeV lower than the lowest spin 1/2 state. Just as we can see in Eq. ( TT7T) . absolute values 
of the non-diagonal matrix elements £°£ E (1,2) and S^S E {2, 1) are more than 2 times of 



those of £°£ E (l,2) and £°S E {2,1), and it's the same case for £f GB (2,3) and £$ GE (3,2), 



'3/2 V x 5 alLKX ^3/2 

1/2 ^ *) miu ^1/2 

so mixing of the spin 3/2 states is much stronger than that of spin 1/2 states as shown in 
Table HVl and W \ and energies of 3/2 states are decreased much more than 1/2 states after 
diagonalization of matrices (|17p . 

In GBE model, the results don't change very much after diagonalization of Eq. (fT7|) . 
Because most of the non-diagonal matrix elements are 0, and the only two nonzero elements 
in Sy BE and £f EE are very small, so mixing of the configurations in GBE model is very weak, 
for instance, for the spin 1/2 states, the only nonzero mixing angle which is for configurations 
2 and 3 is about 0.14. Here we have compared our numerical results to what obtained by 
GBE model in flavor SU(3) symmetric case Jq]. As we can see in Table ITTTj energies of the 
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states with qq in our model are around 100 MeV lower than those obtained by Helminen and 
Riska, and ss are around 220 MeV lower. The reason is that on the one hand the oscillator 
parameter is taken to be ojhr = 228 MeV for all the five-quark configurations in Helminen 
and Riska's model, while in the present manuscript, because there are at least 3 strange 



quarks in the studied configurations, we employ u = ^m/m s UHR — 196 MeV, and take 
into account a correction flH]) caused by lighter constituent mass of the light quark, which 
should contribute about 83 MeV to energies of all the studied configurations, on the other 
hand, considering the corrections of SU (3) symmetry breaking, the strength for a K (or ss) 
exchange between one light and one strange (or two strange) quarks should be smaller than 
strength for ir exchange between two light quarks, but in Ref. the strengths for n, K 
and r\ exchanges are taken to be a same one. 

In INS model, mixing of the configurations with same spin parity quantum number is 
very strong, as shown in Table IIVI and |V] As we can see in Table 1111} energy of the spin 1/2 
state with lowest energy is about 131 MeV lower than the one without configuration mixing 
corrections given in Eq. ( fTTl) . that is caused by strong mixing between first and second 
configurations labeled by Eq. ( Fl5|) and ( Fl6i) . On the other hand, the QQ INS interactions 
lead to mixing between configurations with qq and ss pairs, as shown in Table |V] this kind 
of mixing between the spin 3/2 configurations is very strong. 

Note that it is very difficult to compare our numerical results to the existing experimental 
data, because the spin parity quantum numbers of the experimentally discovered Q* have 
not been clarified. Here we just tentatively compare the present results to what obtained in 
different three-quark models, which are listed in Table [VTl As we can see, the lowest energies 
in the present three five-quark models are ~ 200 MeV lower than most of those listed in 
Table [VTl just as we have expected. Accordingly, five-quark components may be the more 
preferable ones in Q* than traditional three-quark excitations, and this can be examined by 
measurements of Q*Q production in decay performed at BESIII 44J. 



IV. CONCLUSION 

To conclude the present manuscript, here we calculate the energies of 9 low- lying s 3 QQ 
configurations with quantum numbers of Q*. One gluon exchange, Goldstone Boson ex- 
change and instanton-induced interactions are taken as hyperfine interactions between 
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quarks, respectively. In addition, corrections of flavor SU(3) symmetry breaking are also 
taken into account. 

The numerical results show that the OGE and INS hyperfme interactions lead to strong 
mixing of the configurations with same spin parity, and decrease the energies of the studied 
states a lot. In OGE model, it's very interesting that the lowest energy state is with spin 
3/2 but not 1/2, this is different from results obtained by GBE and INS models. In GBE 
model, large discrepancies of our results from what obtained in Ref. |9J indicate that the 
SU(3) breaking corrections are significant. In INS model, QQ hyperfme interactions lead 
to strong mixing between five-quark configurations with light qq and strange ss, this is the 
most significant difference between INS and the other two models. 

It is not convenient for us to compare our results to experimental data, because the data 
is very poor. While comparing to predictions of traditional three-quark models, the state 
with lowest energy in our model is ~ 200 MeV lower, it indicates that the energy cost 
to excite ground state of Q hyperon to a five-quark state is less than that to an orbital 
excitation. Accordingly, five-quark components may be more preferable ones in Q*. If this 
conclusion is correct, namely, the lowest Q* is lying at ~ 1800 MeV, then it can be found 
in measurements of \P' decays performed at BESIII. Moreover, if Q* with lowest energy is 
experimentally found at ~ 1800 MeV, spin parity quantum number of this Q excitation could 
be an evidence to indicate whether the OGE model or the other two hyperfme interaction 
models is more appropriate in studies of Q excitations. 

On the other hand, here we have only considered the configurations with negative parity, 
work on energies of the s 3 QQ configurations with positive parity is in progress. 
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Appendix A: Matrix elements of the operators in OGE model 

Matrix elements of the color operators Af • A?? and — Af ■ A!?* between three color states 

are 

<[211]oi ® CglAf • A^| [211]ci ® Cq) = 4/3 , 
([211] C2 <8 Cg|Af • A£|[211] C2 <g> C ) = -8/3 , 
([211] C3 ® C |Af • A£|[211] C3 <g) C ) = -8/3 , 
([211] C1 <8 Cq| - A^ • Af*|[211]ci (8) Cq) = 2/3 , 
([211] C2 g) C | - A£ • AJ?* |[211] C2 (8) Cq) = 2/3 , 
([211] C3 <g) C | - Af • Xf*|[211]o3 ® C ) = -16/3 . (Al) 

Matrix elements of the spin operators <?i • c? 2 which are independent with the total angular 
momentum of the five quark configurations are 

(S 5 ([22} SI ® Sq)!*! • a 2 \S 5 ([22] sl <g> Sq)) = 1 , 
(£5([22] S2 <g) • a 2 \S 5 ([22] S2 <g> 5 )) = -3 , 

(S 5 ([31] S1 (8 S'q)!^! • <?2|S 5 ([31] 51 ® 5 )) = 1 , 

(S 5 ([31]s2 (8 S-q)!^ ■ ^\S 5 ([31]s2 ® S Q )) = 1 , 
(S 5 ([31] S3 <g) • <? 2 |£ 5 ([31] S3 ® = "3 > 

(S 5 ([4] s ® Sq)!*! • a 2 \S 5 ([A] s ® 5q)) = 1 , (A2) 

and those of <t 4 ■ c? 5 are 

(1/2([31] S1 (8 Sq)|c?4 • cr 5 |l/2([31] si <g) Sq)) = 1 , 
(1/2([31] S2 <8> ^q)|*4 • <? 5 |1/2([31] 52 <g> S Q )) = -2 , 
(l/2([31]s3 ® ^q)I*4 ■ <r 5 |l/2([31] 53 ® 5q)) = -2 , 
(l/2([22] S i ® 5q)|<t 4 ■ <r 5 |l/2([31] 52 (8 Sq)) = \/3, 
(1/2([22] S2 (8 5 )|a 4 • <7 5 |1/2([31] 53 <8> Sq)) = VI , 
(3/2([31] S i (8 5q)|* 4 • <r 5 |3/2([31] 5 i <8> Sq)) = -1/2 , 
(3/2([31] S2 <8 S Q )\a 4 ■ <x 5 |3/2([31] 52 <8> Sq)) = 1 , 
(3/2([31] S3 ® Sq)I*4 • a 5 |3/2([31] 53 ® Sq)) = 1 , 



(3/2([4] 5 <g> S Q )\a A ■ a 5 |3/2([4] 5 ® S Q )) = -3/2 , 
<3/2([31] sl <g> S Q )\a 4 ■ a 5 |3/2([4] s ® 5g)> = >/l5/2 , 

(5/2([4] s <g> S Q )|<7 4 • a 5 |5/2([4] s ® S )> = 1 . (A3) 

Finally, matrix elements of the flavor dependent chromomagnetic strenth operators Ci j2 and 
C 4j 5 are as follow 

(s 3 g([31] F1 ) ® g|C li2 | S 3 g([31] F1 ) <g> g> = |(5C SS + C qs ) , 
(s 3 g([31] F2 ) (g) g|C li2 | S 3 g([31] F2 ) ®q) = \{2C SS + C 9S ) , 

(s 3 g([31] F3 ) ® g|C li2 | S 3 g([31] F3 ) ® ?} = C gs , 
(s 3 g([31] F1 ) (g) g|C lj2 | S 3 g([31] F2 ) ® g> = ^(C 9S - C ss ) , 
( S 3 g([31] F2 ) ® g|C 1)2 | S 3 g([31] F1 ) ® q) = ^{C qs - C ss ) , 
(s 3 q([4} F ) ® g|C li2 | S 3 g([4] F ) ® g) = \{C qs + C ss ) , 
(s\[A] F )®s\C lfl \s\[A] F )®s) = C ss , 
(s 3 g([31] F1 ) ® g|C li2 | S 3 g([4] F ) ® g) = f(C qs - C ss ) , 
(s 3 g([31] F2 ) ® g|C 1 , 2 |s 3 g([4] F ) ® g) = f(C qs - C ss ) , 
(s 3 g([31] F1 ) ® g|C4, 5 |s 3 g([31] F1 ) ® g) = |(3C 3 - g + C gs ) , 
(s 3 g([31] F2 ) ® g|C 4 , 5 |s 3 g([31] F2 ) ® g) = C g - S , 

(s 3 g([31] F3 ) ® g|C4,5|s 3 g([31] F3 ) <8> g) = Q s , 
(s 3 g([4] F ) ® g|C 4 , 5 k 3 g([4] F ) ® g) = ±(C f , + 3Cy , 

(s 4 ([4] F ) ® s|C 4 , 5 | S 4 ([4] F ) ®s) = C- ss ■ (A4) 

Notice that here we have only listed the nonzero matrix elements. 

Appendix B: Matrix elements of the flavor dependent operator in GBE model 

Matrix elements of the flavor dependent operator Cf^Af • A^ follow 

(s 3 g([31] F1 ) ® g|C lj2 |s 3 g([31] F1 ) ® g) = §(10C" S + 2C X ) , 
(s 3 g([31] F2 ) ® g|C lj2 | S 3 g([31] F2 ) ® g) = |(8C SS " + 4C A ') , 

(s 3 g([31] F3 ) ® g|C li2 | S 3 g([31] ra ) ® g) = -§C* 
(s 3 g([31] F1 ) ® g|C li2 | S 3 g([31] F2 ) ® g) = ^(C A - C sS ) , 



(s 3 g([31] F2 ) ® g|C li2 | S 3 g([31] F1 ) ® g> = 2 -f(C K - C*) , 
(s 3 q([A} F ) <g) g|C li2 | S 3 g([4] F ) ® g> = §(C* S + C K ) , 
( S 4 ([4] F )®s|C l!2 | S 4 ([4] F )®s) = |C sS , 
<a 3 g([31] F1 ) ® g|C li2 | S 3 g([4] F ) ® q) = ^{C K - C sS ) , 
(s 3 g([31] F2 ) ® giC li2 | S 3 g([4] F ) ® g> = 4 (C* - C") . (Bl) 



Appendix C: Matrix elements of the operators in the INS model 

Matrix elements of the color operators P^ 6 , P121 P^f anc ^ Pa,5 between three color 
states are 



[211]ci 




pC,6 
^1,2 


[2H]oi 




= 1 , 


[211]c2 




p C,6 
r \,2 


[2H]C2 




= 0, 


[2H]C3 


®C q 


pC,6 


[2H]C3 


®Cq) 


= 0, 


[2H]ci 


®C q 


pC,3 


[2H]ci 


®Cq) 


= 0, 


[2H]C2 


®C q 


p C,3 


[2H]C2 




= 1 , 


[2H]C3 




p C,3 


[2H]C3 




= 1 , 


[2H]ci 


®C q 


p C,8 
-^4,5 


[2H]cx 




= 1 , 


[211] C2 




p C,8 
^4,5 


[211] C2 




= 1 , 


[2H]C3 




p C,8 
■^4,5 


[211]C3 


®Cq) 


= 0, 


[2H]ci 




pC,l 
-"4,5 


[2H]ci 


®Cq) 


= 1 , 


[2H]C2 




pC.l 
-"4,5 


[2H]C2 


®Cq) 


= 0, 


[2H]C3 




|pC,l 
l r 4,5 


[2H]C3 




= 0. 



(CI) 

Matrix elements of the spin operators Pf^ and Pf^ 1 which are independent with the total 
angular momentum of the five quark configurations are 

(S 5 ([22] S1 ® S Q )\Pf^\S 5 ([22]si ® S Q )) = , 
(5 5 ([22] S2 ® 5q)|P 1 5 2 =°|5 5 ([22] 52 ® S Q )} = 1 , 
(S 5 ([22] S1 ® 1 Sq)|P 1 5 2 = 1 |5 5 ([22] 51 ® S Q )) = 1 , 
(S 5 ([22] 52 ® S Q )\P 1 S 2 = 1 \S 5 ([22} S 2 ® S Q )) = , 
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(£ 5 ([31]si$ 




S 5 ([31]si 




= o, 


(5 5 ([31] S2 c> 




S 5 ([31]s2 


®Sq)) 


= 0, 


(ft([31]s3<S 


*S Q )\P§« 


S 5 ([31]s3 


®Sq)) 


= 1, 


(ft([31]si(S 


^Q)IA S 2 =1 


S 5 ([Sl]si 


®Sq)) 


= 1, 


(5 , 5 ([31]52$ 


^Q)IA S 2 =1 


S 5 ([31}s2 


®Sq)} 


= 1, 


(S 5 ([31]s3£ 


^0)IA S 2 =1 


S 5 ([31}s3 


®Sq)} 


= 0, 




^q)IA S 2=° 


S 5 ([4}s® 


Sq)) = 


o, 


(Wise 




\S 5 m s ®S Q )) = 


= 1, 



and those of Pf b and Pf 5 1 are 

(1/2([22] 51 ® ^)|P 4 5 5 = |1/2([22] S1 ® Sq)) = 1/4 , 
(l/2([22] S2 ® S Q )\P 4 s f \l/2([22] S2 ® S Q )) = 1/4, 
(1/2([22] 51 ® 5 )|Pi; 5 = 1 |l/2([22] 51 ® S )> = 3/4 , 
(l/2([22] 52 ® ^)|P 4 s 5 = 1 |l/2([22] 52 ® S^)> = 3/4 , 

(l/2([31] 5 i ® ^)|P 4 5 5 = |l/2([31] 5 i (8) Sq)> = , 
(1/2([31] 52 ® ^)|P 4 5 5 =°|1/2([31] 52 ® fy)) = 3/4 , 
(1/2([31] 53 ® ^)l^4 5 5 =0 |l/2([31] 5 3 ® Sq)) = 3/4 , 

(1/2([31] 51 ® 5g)|Pi; 5 = 1 |l/2([31] 51 ® 5q)> = 1 , 
(l/2([31]s2 ® 5 )|P 4 ? 5 = 1 |1/2([31] 52 (g) S^)) = 1/4 , 

(1/2([31] 53 ® ^)l^4 5 5 =1 |l/2([31]53 ® Sg)> = 1/4 , 
(1/2([22] 51 ® 5 )|Pi; 5 = |l/2([31] 52 ® fy)) = -V3/4 , 
(l/2([22] 52 ® ^)|Pi; 5 = |l/2([31] 5 3 ® Sq)) = -V3/4 , 
(1/2([22] 51 ® 5 )|Pi; 5 = 1 |l/2([31] 52 ® 5g)> = V3/4 , 
(l/2([22] 52 ® 5 )|Pi; 5 = 1 |l/2([31]53® 5 Q )> = v/3/4, 
(3/2([31] 51 ® S )|P 4 5 5 =°|3/2([31] 51 ® Sq)) = 3/8 , 
(3/2([31] 52 ® ^)|P 4 5 5 = |3/2([31] 52 ® = , 
<3/2([31]s3 ® Sg)|P 4 5 5 = |3/2([31]53 ® S Q )) = , 
(3/2([31] 5 i ® 5 )|Pi; 5 = 1 |3/2([31] 51 ® S Q )) = 5/8 , 
(3/2([31] 52 ® 5q)|P 4 5 5 = 1 |3/2([31] 52 ® Sq)) = 1 , 
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<3/2([31] S3 ® ^g)|Pi;5 =1 |3/2([31]53 ® Sq)> = 1 , 
(3/2([4] 5 ® S )|P 4 * 5 =°|3/2([4] 5 s = 5/8 , 

(3/2([4] 5 ® 5 )|Pi; 5 = 1 |3/2([4] 5 ® Sq)} = 3/8 , 

(3/2([31] 51 ®5q)|P 4 5 5 = |3/2([4] 5 ®5q)) = -VTE/8, 

(3/2(i31] sl ®S Q )\Pif\3/2(i4] s ®S Q )) = y/TE/8 , 

(5/2([4] s ® S )|P 4 S 5 =°|5/2([4] S ® S Q )) = , 

(5/2([4] s ® ^)|P 4 5 5 = 1 |5/2([4] 5 ® fy)) = 1 . 
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TABLE I: The studied s 3 QQ configurations with inner quantum number N r = L = 0. Notice that 
Columns FS and CS just denote the same configurations expressed in two different languages. 

FS CS J p 



S 3 q([4] x [21l} c [31} FS [3l} F [22} s ) ® q S 3 g([4] x [31] F [211] cs [211] c [22] 5 ) ® q 

s 3 q([4} x [211] c [31] FS [31] F [31]s) ® q s 3 q([4] x [31] F [211]cs[211] c [31]s) <S) q 

s 3 q([4] x [211]c[3l} FS [31] F [A]s) ® q s 3 q([4] x [31] F [211} C s[2ll]c['i\s) ® q 

s 3 q([4]x[2ll]c[31] FS [4] F [31]s) ® q s 3 q([4] x [4] F [211]cs[211]c[3l]s) ® q 

s 4 ([4] x [211] c [31] F s[4]f[31] 5 ) ® 5 s 4 ([% [4]jr[211] C s[2H]c[31]s) ® 5 
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TABLE II: Values for the parameters in three kinds of hyperfine interactions (in unit of MeV). 



OGE 


111 


340 


771 


460 




196 


Vn 


-208 




Cnr, 


18.3 




11.2 




6.8 




29.8 






18.4 


Css 


8.6 










GBE 


77! 


340 


m s 


460 


U 


196 


V 


-269 




C n 


21 














INS 


m 


340 


m s 


460 


U! 


196 


Vo 


-213 




9 


52.5 


g' 


33.3 











TABLE III: Energies of the studied five-quark configurations in the three models, compared to the 
results obtained in (]] listed in column HR. The results are in unit of MeV. 



J p 


OGE 


GBE 


INS 


HR 


1- 

2 


2146 


1834 


1797 


1917 




2451 


1894 


2027 


1973 




2677 


2012 


2223 


2141 




2655 


2162 


2431 


2381 


3~ 
2 


1818 


1894 


1987 


1973 




2331 


1990 


2025 


2085 




2503 


2012 


2143 


2141 




2517 


2162 


2717 


2381 


5~ 
2 


2492 


1990 


1987 


2085 
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TABLE IV: Coefficients for the mixings of the configurations with spin 1/2 in three models. 



HYP 


1 


2 


3 


4 


OGE 


0.843 


0.538 


0.009 







-0.535 


0.836 


0.122 







0.058 


-0.108 


0.992 
















1 


GBE 


1 
















0.990 


0.138 










-0.138 


0.990 
















1 


INS 


0.678 


0.735 


-0.004 


-0.003 




0.720 


-0.662 


0.184 


0.098 




-0.109 


0.105 


0.949 


0.276 




-0.103 


0.010 


0.256 


0.956 



TABLE V: Coefficients for the mixings of the configurations with spin 3/2 in three models. 



HYP 


1 


2 


3 


4 


OGE 


-0.627 


0.771 


0.110 







0.243 


0.060 


0.968 







0.740 


0.634 


-0.226 
















1 


GBE 


0.990 





0.138 










1 










-0.138 





0.990 
















1 


INS 


0.646 


0.612 


-0.456 


0.000 




0.647 


-0.138 


0.730 


-0.170 




0.352 


-0.560 


-0.254 


0.706 




-0.201 


0.541 


0.440 


0.688 
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TABLE VI: Q* with spin parity 1/2 - and 3/2~ predicted in the non-relativistic quark model 
(OGE), relativized quark model (ROGE), Glozman-Riska model (GBE), covariant quark model 
based on Bethe-Salpeter equation with instanton-induced quark force (BSE), large N c analysis, 
algebraic model (BIL), and Skyrme model. The mass is given in the unit of MeV. 

State OGE [6j ROGE [7j GBE [8j BSE[10j Large-iV c [11-.15J BIL [43] Skyrme Model[16j 
2020 1950 1991 1992 2061 1989 1837 

2410 2456 
2020 2000 1991 1976 2100 1989 1978 

2440 2446 2604 
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